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ABSTRACT 

We propose to use relative galaxy ages as a means of constraining cosmological pa- 
rameters. By measuring the age difference between two ensembles of passively-evolving 
' galaxies at somewhat different redshifts, one could determine the derivative of redshift 

■ with respect to cosmic time, dz/dt. At high redshifts, z ~ 1-2, this measurement would 

^ . constrain the equation-of-state of the dark energy, while at low redshifts, z < 0.2, it 

^ . would determine the Hubble constant, Hq. The selected galaxies need to be passively- 

, evolving on a time much longer than their age difference. 
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Subject headings: Cosmology: theory — galaxies 



Introduction 



Recently, there has been much progress in constraining cosmological parameters, such as the 
matter content of the Universe (Peacock et al. 2001) and the Hubble constant (Preedman et al. 
2001). Available data for the microwave background anisotropics on degree scales (de Bernardis et 
^ , al. 2000; Hanany et al. 2000; Netterfield et al. 2001; Lee et al. 2001 ) and the Hubble diagram 

Y^. ' of Type la supernovae (Riess et al. 1998; Perlmutter et al. 1999) indicate that the Universe 

^ ■ has a flat geometry and is dominated at present by some form of dark energy with a negative 

pressure (Garnavich et al. 1998; Perlmutter et al. 1999). The equation of state of the dark energy, 
Pq = wqpq, expresses the ratio between the pressure, pq, and the mass density, pg, of the dark 
energy in terms of the parameter wq (in units of c = 1). The value of wq could either be constant as 
I in the case of a cosmological constant {wq = —1), or time-dependent as in the case of a rolling scalar 

field or "Quintessence" (Ratra & Peebles 1988; Caldwell et al. 1998). Any such behavior would 
have far-reaching implications for particle physics. Therefore, the next observational challenge is 
to determine the evolution of wq as a function of redshift (Huterer & Turner 2000; Maor et al. 
2001; Weller &; Albrecht 2001). The related observations need to be done at redshifts z <2, when 
the dark energy started to dominate the expansion of the Universe. 

The popular approach for measuring wq{z) uses its effect on the luminosity distance of sources. 
In particular, the proposal for the Supernova/ Acceleration Probe (SNAP) mission^ suggests to 
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monitor < 2000 Type la supernovae across the sky per year and determine their luminosity distances 
up to a redshift z ^ 1.5 with high precision. However, the sensitivity of the luminosity distance to 
the redshift history of wq{z) is compromised by its integral nature (Maor et al. 2001), 



where t{z) is the age of the Universe at a redshift z which depends on wq{z). 

In this paper, we propose an alternative method that offers a much better sensitivity to wq{z) 
since it measures the integrand of equation (1) directly. Any such method must rely on a clock 
that dates the variation in the age of the Universe with redshift. The clock in our method is 
provided by spectroscopic dating of galaxy ages. Based on measurements of the age difference, At, 
between two passively-evolving galaxies that formed at the same time but are separated by a small 
redshift interval Az, one can infer the derivative, (dz/dt), from the ratio (Az/At). The statistical 
significance of the measurement can be improved by selecting fair samples of passively-evolving 
galaxies at the two redshifts and by comparing the upper cut-off in their age distributions. All 
selected galaxies need to have similar metallicities and low star formation rates (i.e. a red color), 
so that the average age of their stars would far exceed the age difference between the two galaxy 
samples. At. 

This differential age method is much more reliable than a method based on an absolute age 
determination for galaxies (e.g., Dunlop ct al. 1996; Alcaniz & Lima 2001; Stockton 2001). As 
demonstrated in the case of globular clusters, absolute stellar ages are more vulnerable to systematic 
uncertainties than relative ages (Stetson, Vandenberg &; Bolte 1996). Moreover, absolute galaxy 
ages can only provide a lower limit to the age of the Universe and only place weak constraints on 
the possible histories of wq{z). 

The quantity measured in our method is directly related to the Hubble parameter. 



Hence, an application of this method to old elliptical galaxies in the local Universe can be used to 
determine the value of the current Hubble constant, Hq. 

In §2 we quantify the advantage of the above differential age method relative to the standard 
luminosity distance method in constraining the history of wq{z). We then describe the spectro- 
scopic dating technique and apply it to mock galaxy spectra in §3. In particular, we analyze the 
dependence of the cosmological constraints on the signal-to-noise ratio of the spectroscopic data 
and the number of observed galaxies. The significance level of the constraints attainable with a 
single pair of galaxies, dictates the number of such pairs required in order to differentiate between 
various histories of wq{z). Finally, we summarize our main conclusions in §4. 
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Observables 



We consider a flat universe composed of matter and dark energy, the latter having an equation 
of state pQ = wq{z)pq, where wq may depend on redshift. The Hubble parameter is = 
Hq[pt{z)/pt{0)]. Here, the subscripts Q, m, and T refer to the dark energy, the matter, or 
the total sum of the two, respectively. Assuming further that the matter is non-relativistic (i.e. 
effectively presureless), we get 
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where we have used the energy conservation equation for the dark energy, pq = — 3i?(l + wq)pq 
(Maor et al. 2000). Thus, (dz/dt) is related to the equation of state of the dark energy through 
one integration only, while the luminosity distance in equation (1) is given by an integral of the 
inverse of equation (3), namely through two integrations. By differentiating (dz/dt) with respect 
to t we find, 
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which depends explicitly on wq without any integrations. Thus, the second derivative of redshift 
with respect to cosmic time measures wq directly. While it is possible to find significantly different 
redshift histories of wq for which the evolution of o?l is similar, this cannot be done for [d^z/dt^). 

Maor et al. (2000) have argued that due to its integral form, the luminosity distance has only 
a weak discriminating power with respect to different possible histories of wq{z). To demonstrate 
their case, they considered examples of wq which have a quadratic dependence on z. Figure 1 shows 
similar examples (top panel) along with their corresponding observables, namely the luminosity 
distance di,{z) (second panel), dz/dt (third panel), and (Pz/dt^ (bottom panel). The different 
wq{z) histories generate a variation of ~ 5% in t^L, ~ 10% in (dz/dt), and ~ 30% in ((Pz/dt^) at 
z ~ 1.5. The last observable has a fractional variation as large as that of wq(z). 



3. Spectroscopic Dating of Passively— Evolving GallELxies 

The change of cosmic time with redshift may be inferred from the aging of stellar populations 
in galaxies. This inference must be done with caution since galaxies are vigorous sites of star 
formation at high redshift. It might seem difficult to estimate accurately the differential aging of 
the universe based on star-forming galaxies, since the stars in these galaxies are born continuously 
and a young stellar population may dominate their emission spectrum. Fortunately, examples of 
passively-evolving galaxies have already been identified in large numbers at high redshifts (Dunlop 
1996; Daddi, Cimatti & Renzini 2000; Dey et al. 2001; Stockton 2001). We focus our attention on 
these red galaxies since their light is dominated by an old stellar population. 
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Our goal is to measure the age difference between two passively-evolving galaxies at different 
redshifts. How accurately can this be done? We emphasize that the relative age is better determined 
than the absolute age since systematic effects on the absolute scale are factored out (for a frac- 
tional age difference <C 1). At redshifts z ~ 1-2 the rest- frame UV spectrum of non star-forming 
(elliptical) galaxies is dominated by light from main sequence stars with masses of 1-2 Mq, in a 
regime of stellar evolution that is well understood (Spinrad et al. 1997). The lack of interstellar 
dust or gas simplifies further the spectrum of old elliptical galaxies. The top panel in Figure 2 
shows a simulated galaxy spectrum in the rest frame UV that can be easily collected by 10 meter 
class telescopes from the ground (e.g. Dunlop ct al. 1996). This spectrum was calculated for a 
single stellar population with a specific metallicity (solar) and age (2.4 Gyr), using the synthetic 
stellar population models developed in Jimenez et al. (1998). The simulated spectrum has a 10 A 
resolution, to which we have added Poisson photon noise with a signal-to-noise ratio S/N = 30 per 
resolution element. The bottom panel shows the values recovered at the 95% confidence level (using 
a goodness-of-fit statistics) for the metallicity and age with S/N = 30 (dark shaded region) and 
S/N = 50 (lightly shaded inner region). Note that age and metallicity are not degenerate in the 
rest-frame UV, provided that the measured spectra have a sufficient S/N (Nolan et al. 2001). A 
value of S/N = 15-30 can be obtained for a faint galaxy of magnitude 24 in the i?-band after 10 
hours of integration on the Keck telescope (K. Adelberger, private communication). Observation of 
a sufficiently wide field of view would allow to measure spectra of many passively-evolving galaxies 
at the same time. 

The assumption that stars in an elliptical galaxy have a single metallicity value is clearly 
an over-simplification. Local ellipticals have metallicity gradients (e.g. Davies et al. 1993), even 
though these gradients are mild (e.g. Friaca & Terlevich 1998). We therefore consider the possibility 
that elliptical galaxies have a mixture of populations with different metallicities*^ . In order to assess 
the accuracy by which we can recover the age of a mixed stellar population, we simulated several 
spectra with a mixture of six metallicity values ranging from 0.01 to 5 times the solar value. For 
example, one of our models had weights of 5, 10, 25, 20, 12 and 28% for metallicities of 0.01, 0.2, 1, 
2.5 and 5 times the solar value, respectively. We then attempted to recover the age by searching for 
the best-fit model with a single metallicity, and found an age uncertainty of 0.1 Gyr. Alternatively, 
we fitted the simulated spectrum by a mixed metallicity model with six free components, and 
recovered a more accurate age with an uncertainty of ~ 0.05 Gyr (as well as a reduced metallicity 
uncertainty of 8%). 

Using the more conservative age uncertainty of 0.1 Gyr derived above for S/AT = 30, we may 
now estimate the attainable constraints on wq{z) for a fiducial pair of passively-evolving galaxies 
which formed at the same time but are observed at 2; = 1.4 and 1.6 and have ages of 2.79 and 



*The small spread in age among stars in an elliptical galaxy does not affect our results, since the method only 
relies on the shift in the average stellar age as a function of redshift. For a passively-evolving galaxy, the shift in the 
average age of its stars within a given redshift interval, is determined only by the background cosmology. 
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2.4 Gyr, respectively. Figure 3 shows the resulting error bar on (dz/dt), which is not sufficiently 
restrictive for a single pair of galaxies. However, for a statistical ensemble of independent galaxies 
(which do not all belong to the same local region of a galaxy cluster, for example), the error bar 
will be reduced by the square root of the number of analyzed galaxy pairs. The smaller error bar in 
the plot illustrates the constraint that could be placed by 20 such pairs. This constraint is sufficient 
to distinguish among the different wq{z) histories in Figure 1. 

The situation in reality is not as simple as described above, since a galaxy survey will not 
just identify galaxies that formed at the same cosmic time. Rather, the galaxies selected as not 
having substantial star formation, will show an age distribution at any given redshift. Given this 
age distribution, what is the accuracy by which one may recover the age shift, provided that the age 
of each passively-evolving galaxy is measured with a 0.1 Gyr uncertainty? To answer this question, 
we have simulated the galaxy age distributions at two different redshifts. The distribution at the 
higher redshift was modeled fiducially as a Schechter function which is asymptotically flat at low 
ages and is fully truncated at some finite maximum age value (because galaxies cannot be older 
than the age of the universe). In our example we chose a cutoff at 2.4 Gyr. The distribution 
was then binned into bins of 0.1 Gyr width, corresponding to the age uncertainty for individual 
galaxies. The age distribution at a lower redshift was chosen to be the same as the first one 
but shifted by either 0.0, 0.15 or 0.30 Gyr (so that the maximum age in the distribution would 
be 2.4, 2.55 or 2.7 Gyr, respectively). We also removed or added a fixed fraction (20%) of the 
galaxies within each bin through a random Poisson realization of this fraction. This degree of 
freedom is required since some galaxies may have merged with a star-forming galaxy and therefore 
disappeared from the old galaxy population. Moreover, the normalization of the two distributions 
is regarded as a free parameter because the volume surveyed at the two redshifts might be different. 
Our approach assumes that both samples of galaxies have been passively evolving for a time much 
longer than the age difference between the two redshifts. Each member of the low-redshift sample 
has a progenitor which, statistically speaking, is included in the high-redshift sample. But each 
high-redshift galaxy has some probability of encountering a merger and disappearing from the 
low-redshift sample. Our method would fail only if this probability is both significantly large 
and significantly age-dependent. However, for a sufficiently narrow redshift interval, the merger 
probability is small. The tight color-magnitude relation of present-day ellipticals further limits the 
effect of mergers on their age distribution (Peacock 1991; Bower, Lucey Sz Ellis 1992; Peacock et 
al. 1998). 

The two distributions are shown as the two upper panels of Figure 4 for an age shift of 0.3 
Gyr and a total of ~ 75 galaxies in each distribution. We binned the distributions for illustrative 
purposes. The best statistical method to recover the shift is the Kolniogorov-Smirnov (KS) test 
on the unbinned data. The bottom panel shows the probability for getting a particular value of 
the peak in the KS probability (Press et al. 1992) for the three cases wc consider. Wc derived 
this probability distribution through a large number of Monte-Carlo realizations of the above 
distributions. Also shown as shaded regions are the shift levels for which one could distinguish 
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among the different wq{z) histories in Figure 3. We find that the shifts can be recovered with good 

precision, despite the fact that the tip of the age distribution is populated by a small number of 
galaxies. In particular, we get a > 90% probability in this example for distinguishing among the 
different wq{z) histories. 

At present, there is no such set of galaxy spectra that could discriminate between a constant 
and a variable wq{z). Nevertheless, it is instructive to examine the best spectroscopically observed 
pair of red galaxies to date, namely 53W091 and 53W069, which reside at z = 1.43 and 1.55, 
respectively. The stellar populations of these galaxies have been studied in detail recently by Nolan 
et al. (2001), who demonstrated that the age-metallicity degeneracy can be broken based on UV 
rest-frame spectroscopy with high quality. However, since these galaxies where observed with only 
(S/N) ~ 6-10, they do not provide strong constraints on wq{z). 

4. Conclusions 

We have shown that the derivative of cosmic time with respect to redshift can be used as an 
effective tool for measuring the equation of state of the dark energy as a function of redshift. For 
this purpose, spectra of passively-evolving galaxies need to be obtained with a high signal-to-noise 
ratio. In order to differentiate at the 2-a level between a constant and a variable wq{z), one needs 
~ 70 pairs of galaxies at i? ~ 24. For a field of view ^ (10')^, which is attainable with instruments 
like DEIMOS on the Keck telescope and VIMOS on the VLT, we estimate that the required number 
of ellipticals may be found in the volume between z = 1.4 and 1.6, and so all pairs can be observed 
simultaneously over 10 hours of integration on a 10 meter class telescope, such as Keck or the VLT. 

The key advantages of our method relative to supernova searches such as the SNAP space 
mission, are that (i) the observable (dz/dt) is more sensitive to wq{z) than the luminosity distance 
c^l; and (ii) the related precision observations can be done from the ground. In fact, some observa- 
tional groups have already identified relevant sets of red galaxies at high-redshifts (e.g., Daddi et 
al. 2000; Stockton 2001). 

The Hubble constant, Hq, can be measured by computing the age difference between red 
galaxies at z = and z ~ 0.2. The characteristic age difference in this redshift range, ~ 2 Gyr, can 
be measured for a single pair of galaxies with a precision of 20% given data with S/N = 10 and 
spectral resolution of lOA . A large volume of spectroscopic data with this quality on red galaxies in 
the local universe, will soon be released by the Sloan Digital Sky Survey (SDSS^). The wavelength 
coverage of the SDSS spectra of 3900-9100 A includes spectroscopic features such as il^, that are 
sensitive to the aging of the stellar population. Given a statistical sample of many hundreds of red 
SDSS galaxies (D. Eisenstein, private communication), one could in principle determine the value 
of Ho to a percent accuracy, as long as the systematic modeling uncertainties for relative ages can 



^ http : / / www. sdss . org 
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be reduced to that level. The differential age method is independent of the Cepheid distance scale 
and is subject to other uncertainties. 

We thank Kurt Adclbcrgcr, Rcnnan Barkana, Rob Caldwell and Marc Kamionkowski for useful 
discussions. This work was supported in part by NASA grants NAG 5-7039, 5-7768, and by NSF 
grants AST-9900877, AST-0071019. 
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Fig. 1. — Several examples of wq{z) as a function of redshift (top panel) and the sensitivity of 
several observables: c/l (second panel), dz/dt (third panel) and (Pz/dt^ (bottom panel). 
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Fig. 2. — The top panel shows a simulated spectrum (flux per unit wavelength, Fx) for a galaxy 
with a solar metallicity and an age of 2.4 Gyr, including Poisson photon noise with S/N = 30 
at a spectral resolution of lOA. The bottom panel shows the accuracy of the metallicity and age 
recovery for S/N = 30 (dark shaded region) and S/N = 50 (lightly shaded region). 
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Fig. 3. — The large error bar illustrates the constraint that can be placed on wq{z) by using only 
a single pair of red galaxies that formed at the same time and are observed at z = 1.6 and z = 1.4. 
The spectra are assumed to have a spectral resolution of lOA and S/N = 30. When data for 20 
such pairs is available, the uncertainty is reduced to the thicker, smaller error bar. 
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Fig. 4. — The top panel shows a hypothetical age distribution for 75 galaxies at a redshift z = 1.6. 
The middle panel shows the same distribution at z = 1.4, in which all galaxies are older by 0.3 
Gyr and their fraction per bin changes by ±20% with Poisson fluctuations (due to mergers). The 
bottom panel shows the probability distribution of the recovered shift in age between the two galaxy 
distributions using the peak probability value of the KS-likelihood test for unbinned data with an 
age error of 0.1 Gyr per galaxy. We show three cases of shifts, namely: 0.0, 0.15 and 0.3 Gyr. The 
shaded regions show the range in the recovered shift which would allow to separate the different 
wq{z) histories in Figure 3. 



